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Numerical results from large N reduced QCD2 
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Some results in QCD2 at large N are presented using the reduced model on the lattice. Overlap fermions are 
used to compute meson propagators. 



1. Introduction 

This follows the talk given by H. Neuberger [Q 
and summarizes the numerical results obtained 
in QCD 2 at large N using the quenched reduced 
model and overlap fermions 

The reduced model has a single site, and there 
are d SU(N) matrices U^, fi = 1, • • • , d. The gauge 
action is given by 



s q = -p 



= \U„U V ], (1) 



where (3 = ^p- and g 2 is the continuum gauge 

coupling. We will keep b 2 — ^ fixed as N goes to 
infinity, and this amounts to setting g 2 N = ^p-. 

The reduced model on the lattice has a [[/(l)]" 1 
symmetry. The Eguchi-Kawai argument for re- 
duction holds only if the symmetry re- 
mains unbroken ||. This is not the case as 
b 2 — > 00 H . This problem is resolved by defining 
the quenched reduced model Q where the eigen- 
values of Un, 



a -0-5 



diag(e l 



),(2) 



are fixed, and V^s are the only dynamical degrees 
of freedom. The 9^ are randomly distributed in 
the interval [— n, tt], and a quenched average over 
these variables is taken. However in QCD2, it is 
sufficient to pick 8 l to be the N roots of unity, 
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Figure 1. Plot of P L {Q 2 ,l) at N c = 41 for 
b 2 = 1,2,3,4,5. 
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Figure 2. Plot of the meson propagator for mass- 
less quarks as a function of b 2 . The b 2 = 1 data 
is at N c = 37, the b 2 = 2, 4 data are at iV c = 41 
and the b 2 = 3, 5 data are at AT r = 47. 



and there is no need to perform a quenched av- 
erage Q. Minima of the action occur when all 
the U^s are simultaneously diagonal. Certain 
V^s relate one minimum with another minimum. 
Numerical simulations should sample all minima 
properly. 

Fermionic loops are suppressed by 1/N in the 
large N limit, and only the gauge degrees of free- 
dom determine the dynamics. One can compute 
the fermionic propagator in momentum space by 
force feeding momenta [^) . This amounts to using 
e lp i* Up for the fermion parallel transporter in the 
fi direction which carries a momentum p^. Wc 
will use the propagator for massive overlap Dirac 



Figure 3. Plot of the current-current correlator 
at 7 = 1 as a function of b 2 at N = 41. 



operator given by |6| 

1 1 - jd+ie(H w ) 



g(p) 



2m (1 + fi) + (1 - (j,)~/ d+1 e(H w ) ' 



(3) 



The quark mass is m q — 2m/i where m is the mag- 
nitude of the negative Wilson mass that appears 
inside H w , the hermitian Wilson-Dirac operator. 

2. Correlation functions 

Meson correlation functions are known analyt- 
ically in QCD2 0], and we can check the results 
obtained from the quenched reduced model on 
the lattice. Correlators can be expressed as an 
infinite sum over stable mesons, which reproduce 
free field behavior at large momenta. Therefore 
we have confinement and asymptotic freedom in 



3 



QCD2. For m q — 0, there is a massless meson 
in the spectrum, which is interpreted as a Gold- 
stone boson. Chiral symmetry breaking occurs if 
we take m q — » after we take N — > oo. 

Gauge fields are generated using (||) and the 
quenching condition in (^) . The eigenvalues of U\ 
and U 2 are set to the N roots of unity, and gauge 
field configurations are generated at a fixed b 2 for 
different values of N. Fermionic momenta also are 
fixed at N roots of unity. The "bare" scalar and 
pseudoscalar propagators are computed using 



So(P) 



— y 

AN 2 ^ 

P1.P2 

Tr[g(pi + Pi,p 2 + Pb)g(pi,P2) 

g{pi - Pi,p 2 + P 2 )g(pi,p 2 ) 
g{pi + P2,P2 + Pi)g(pi,p 2 ) 
g{pi + P 2 ,p 2 - Pi)g(pi,p 2 )} 



Po(P) 



1 



E 



AN 2 

P1,P2 

Tr[g(p 1 + P 1 ,p 2 + P 2 )g^p 1 ,p 2 ) 

+ gipi - Pi,p 2 + P 2 )gHpi,p 2 ) 
+ g(pi + P 2 ,p 2 + P\)g\pi,p 2 ) 
+ g(p\ + P 2 ,p 2 - P\)g\pi,p 2 )] 

This results in a correlation function for mesons 
that is a function of one variable, namely, P 2 = 
4 (sin 2 -EL+sin 2 -&■). The scale is set by e 2 = 2^2, 

and we will use 7 = as a dimensionless mea- 
sure of the quark mass. We will also use a di- 
mensionless momentum Q 2 = The bare cor- 
relators have to be regularized. This is done 
by a subtraction at zero momentum for massive 
quarks and a subtraction at non-zero momentum 
for massless quarks. We will refer to the renor- 
malized quantities on the lattice by Sl(Q 2 ,j) 
and Pl(Q 2 , 7). In the dimensionless notation, we 
have to compare these results on the lattice with 
the regularized results in the continuum 



P R (Q\l)= £ 



n>0 even 



> 2 + /4 



n>l odd 



Q 2 



/'n 



(4) 



(5) 



The dimensionless meson masses fj, n are obtained 
by solving the 't Hooft Hamiltonian ^ , and /i 2 ~ 
7r 2 n in the asymptotic limit. The residues, r 2 , are 
given by 



tlx 



(6) 



where <fi n (x) are the eigenvectors of the 't Hooft 
Hamiltonian. 

Figure |l| shows the results for Pl (Q 2 ,1) for var- 
ious values of b 2 and the comparison to the contin- 
uum result. One can also compute the correlation 
function for massless quarks at finite N. In this 
case, one will not see the effect of chiral symme- 
try breaking. The propagators for the scalar and 
the pseudoscalar are the same and should match 
with the average of the two in the continuum. 
Figure ^| shows the results for massless quarks 
along with the continuum result. One can also 
compute current-current correlators and obtain 
n(Q 2 ). We present the results for current-current 
correlators at 7 = I along with the continuum re- 
sult in figure |^. 

3. Conclusions 

The numerical results obtained in QCD 2 using 
the quenched reduced model and overlap fermions 
are promising. The next step is to perform the 
same computation for QCD4, and this project is 
under way. 
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